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ABSTRACT 

The A cold dark matter (ACDM) standard cosmological model is in severe tension with several cosmological observations. 
Foremost is the Hubble tension, which exceeds 50 confidence. Galaxy number counts show the Keenan—Barger—Cowie (KBC) 
supervoid, a significant underdensity out to 300 Mpc that cannot be reconciled with ACDM cosmology. Haslbauer et al. 
previously showed that a high local Hubble constant arises naturally due to gravitationally driven outflows from the observed 
KBC supervoid. The main prediction of this model is that peculiar velocities are typically much larger than expected in the 
ACDM framework. This agrees with the recent discovery by Watkins et al. that galaxies in the CosmicFlows-4 catalogue have 
significantly faster bulk flows than expected in the ACDM model on scales of 100 — 250 h~! Mpc. The rising bulk flow curve 
is unexpected in standard cosmology, causing 4.80 tension at 200h~! Mpc. In this work, we determine what the semi-analytic 
void model of Haslbauer et al. predicts for the bulk flows on these scales. We find qualitative agreement with the observations, 
especially if our vantage point is chosen to match the observed bulk flow on a scale of 50h~! Mpc. This represents a highly 
non-trivial success of a previously published model that was not constrained by bulk flow measurements, but which was shown 
to solve the Hubble tension and explain the KBC void consistently with the peculiar velocity of the Local Group. Our results 
suggest that several cosmological tensions can be simultaneously resolved if structure grows more efficiently than in the ACDM 
paradigm on scales of tens to hundreds of Mpc. 
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1 INTRODUCTION 


The Hubble tension is one of the greatest currently debated unsolved 
problems in cosmology (Wong et al. 2020; Migkas et al. 2021; Di 
Valentino 2022; Riess et al. 2022). It is a statistically significant 
discrepancy between direct local measurements of the Hubble- 
Lemaitre constant Ho and the prediction of the Lambda-Cold Dark 
Matter (ACDM; Efstathiou, Sutherland & Maddox 1990; Ostriker & 
Steinhardt 1995) standard model of cosmology with parameters 
calibrated to fit the angular power spectrum of anisotropies in the 
cosmic microwave background (CMB). The local Universe appears 
to be expanding 10 percent faster than this prediction. While the 
origin of this tension is not known, if the high local determination of 
Ho is correct, then the universe would have to be about 10 per cent 
younger than if the lower CMB-based value is correct. However, the 
ages of the oldest stars argue against this interpretation (Cimatti & 
Moresco 2023). While their upper limit on Ho is consistent with CMB 
measurements taken by the Planck satellite (Planck Collaboration VI 
2020) and the Atacama Cosmology Telescope (Aiola et al. 2020), 
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it leaves open the issue of why so many local measurements show 
a faster expansion rate (e.g. figure 1 of Di Valentino 2021, and 
references therein). It has recently been argued that the Hubble 
tension should be solved largely at late times in cosmic history 
(Jia, Hu & Wang 2023; Vagnozzi 2023), with the expansion rate 
apparently diverging from ACDM expectations only rather recently 
(Gémez-Valent et al. 2023). 

Another serious but less widely known problem is the Keenan— 
Barger-Cowie (KBC) void, an underdensity with a diameter of 
about a Gpc (Keenan, Barger & Cowie 2013). Such an extended 
and deep local underdensity is indicated by several observations at 
optical (Maddox et al. 1990), infrared (Huang et al. 1997; Frith 
et al. 2003; Busswell et al. 2004; Frith, Shanks & Outram 2005; 
Frith, Metcalfe & Shanks 2006; Keenan, Barger & Cowie 2013; 
Whitbourn & Shanks 2014; Wong et al. 2022), X-ray (Bohringer 
et al. 2015; Böhringer, Chon & Collins 2020), and radio wavelengths 
(Rubart & Schwarz 2013; Rubart, Bacon & Schwarz 2014). The near- 
infrared measurements imply that the matter density is only half the 
cosmic mean value out to a distance of 300 Mpc (see figure 11 of 
Keenan, Barger & Cowie 2013 and fig. 1 of Kroupa 2015). Using 
data from the Millennium XXL simulation (MXXL; Angulo et al. 
2012), Haslbauer, Banik & Kroupa (2020, hereafter HBK20) showed 
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that the KBC void is in 6.040 tension with the ACDM model despite 
accounting for redshift space distortions induced by outflow from 
the void, which implies that the actual density contrast is about half 
that reported by Keenan, Barger & Cowie (2013). Such a deep and 
extended void suggests a cosmological model in which structure 
grows faster than in ACDM. 

It has been suggested that outflows from a large local void can 
solve the Hubble tension (Keenan, Barger & Cowie 2016; Shanks 
et al. 2019; Ding, Nakama & Wang 2020; Camarena et al. 2022; San 
Martin & Rubio 2023). This would be a natural solution because 
a local void creates a hill in the potential, causing galaxies to flow 
away from the void. From a purely kinematic perspective, the near 
uniformity of the CMB implies that a substantial underdensity today 
must be a consequence of outflows. Using this argument, HBK20 
estimated that the observed KBC void implies that the locally 
measured Hy should exceed the true background expansion rate 
by 11 percent (see their section 1.1). The Local Group (LG) is 
situated within the KBC void and is a part of the bulk matter flow. 
While the Hubble tension can be naturally understood as arising from 
outflows induced by the KBC void, the problem for ACDM is that 
it is unable to form such a void because on the relevant scale of 
300 Mpc, the universe should be almost homogeneous and isotropic. 
A simultaneous solution to the KBC void and Hubble tension is 
possible in ACDM, but only if we assume a 100 density fluctuation 
at early times (see figure 1 of HBK20). Moreover, the observations 
of Keenan, Barger & Cowie (2013) trace the majority of the galaxy 
luminosity function, suggesting that the observed underdensity is 
a genuine underdensity in the total matter distribution (see their 
figure 9). Allowing for redshift space distortions, HBK20 found that 
the results of Keenan, Barger & Cowie (2013) imply that the density 
within 300 Mpc is ~ 20 per cent below the cosmic mean evident in 
galaxy number counts at greater distances. This is very much in line 
with the findings of Wong et al. (2022). 

HBK20 showed that matter inhomogeneities comparable to the 
KBC void can arise in the neutrino Hot Dark Matter (vVHDM) 
cosmological framework (Angus 2009; Katz et al. 2013; Wittenburg 
et al. 2023). The vHDM model assumes Milgromian dynamics 
(MOND; Milgrom 1983, 2014). MOND postulates that the gravi- 
tational acceleration g in an isolated spherically symmetric system 
is asymptotically related to the Newtonian gravitational acceleration 
8y Of the baryons alone according to 
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The key new ingredient is a fundamental acceleration scale a, = 
1.2 x 107!° m s~?. This value must be deduced empirically, just 
like the gravitational constant G in standard gravity. Different 
studies over the decades have returned very similar values for a, 
(Begeman, Broeils & Sanders 1991; Gentile, Famaey & de Blok 
2011; McGaugh, Lelli & Schombert 2016). The MOND gravitational 
field follows from a Lagrangian, ensuring the usual symmetries and 
conservation laws with respect to the linear and angular momentum 
and the energy (Bekenstein & Milgrom 1984; Milgrom 2010). The 
least action principle then leads to a generalized Poisson equation that 
is non-linear in the mass distribution. 

MOND has been highly successful on galaxy scales (Famaey & 
McGaugh 2012; McGaugh 2020; Banik & Zhao 2022). vHDM 
extends its application to cosmology by postulating an extra species 
of sterile neutrino with a rest energy of 11 eV, which is crucial to 
matching the high velocity dispersions of virialized galaxy clusters, 
the offset between the weak lensing and X-ray peaks in the Bullet 
Cluster, and the high third peak in the angular power spectrum of the 
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CMB (for a review, see section 9.2 of Banik & Zhao 2022). The back- 
ground cosmology in vHDM is the same as in ACDM (Skordis & 
Ziosnik 2019), though the extra neutrino species implies a mild 
departure from standard Big Bang nucleosynthesis (see section 3.1.2 
of HBK20). Moreover, both models behave similarly at early times 
because when the redshift z = 50, the typical accelerations are high 
and thus standard gravity applies. At lower redshifts, significant 
differences arise in the large-scale structure due to the lack of 
CDM and the MOND corrections to the gravitational field (Angus & 
Diaferio 2011; Angus et al. 2013; Katz et al. 2013; Wittenburg et al. 
2023). 

Using a semi-analytic model in the vHDM framework, HBK20 
evolved three initial void density profiles (Maxwell—Boltzmann, 
Gaussian, and Exponential) for a large grid of initial void sizes 
and strengths at z = 9. The initial conditions were constrained 
by observations of the local Universe (z = 0.01-0.15) and the 
requirement for the density to almost recover the cosmic mean value 
at distances of 600-800 Mpc (see fig. 11 and table 1 of Keenan, 
Barger & Cowie 2013). Any local void solution to the Hubble tension 
implies quite large peculiar (CMB-frame) velocities, so an important 
constraint on such models is the observed peculiar velocity of the LG 
(vto = 627 + 22 km s™!; Kogut et al. 1993). While peculiar velocities 
are typically larger in VHDM (Katz et al. 2013; Wittenburg et al. 
2023), such a low value arises for a reasonable fraction of observers 
because gravitational fields from nearby and more distant structures 
can sometimes partially cancel. However, if we consider a larger 
region of the universe, the significant Milgromian enhancement to 
gravity implies much more substantial bulk flows of galaxies on 
scales of hundreds of Mpc. If this model is to be viable, such rapid 
bulk flows need to be verified observationally. 

We test the local void solution to the Hubble tension proposed 
in HBK20 by extracting the predicted bulk flow from the exact 
same model without any adjustments in order to compare its a 
priori bulk flow predictions with recently published measurements. 
Using the CosmicFlows-4 galaxy catalogue (Tully et al. 2023), 
Watkins et al. (2023) present the bulk flow of galaxies on scales 
of 100 — 250 h~! Mpc, where h ~ 0.7 is the Hubble constant in units 
of 100 km s~! Mpc™!, velocities are reported in the CMB frame, 
and the bulk flow involves a vector average of line of sight peculiar 
velocities out to some distance (Section 2.3). Bulk flows are expected 
to decrease with scale, but the observed bulk flow curve has the 
opposite behaviour and rises to >400 km s7! beyond 16047! Mpc. 
At a scale of 200 h7! Mpc, it is in 4.810 tension with the ACDM 
model (P = 1.49 x 10~°).! An independent study recently reported 
‘excellent agreement with the bulk flow measurements of Watkins 
et al. (2023)’ using the same data set and also found significant 
tension with ACDM at an effective depth of 173h~! Mpc, though 
the more conservative methodology prevented the authors from 
going further out (Whitford, Howlett & Davis 2023). This issue 
is observationally unrelated to the Hubble tension because adopting 
a different Hp would affect the peculiar velocities in a spherically 
symmetric manner, thus not affecting the inferred bulk flow curve 
shown in fig. 7 of Watkins et al. (2023). In Section 3, we compare 
the results in the bottom right panel of this figure with the bulk flow 
predicted by the vHDM model for different void density profiles 
and possible vantage points that were previously shown to provide 
the best overall match to several other cosmological observables 


'The tension is not stated for a scale of 250h~! Mpc, but figure 8 of 
Watkins et al. (2023) shows that it exceeds the 5o falsification threshold 
(P = 5.73 x 1077) typically used in science. 
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(mainly the KBC void density profile and the high local Hp; 
see HBK20). 

It is important to note that the bulk flow measurements of Watkins 
et al. (2023) were not available to HBK20. While earlier bulk flow 
measurements were available, these were not used — the vHDM 
model and the parameters of a local void in this model were in no way 
constrained to fit the observed bulk flow on scales of several hundred 
Mpc. The primary objective of this study is to use the bulk flow 
measurements of Watkins et al. (2023) to test the a priori predictions 
of the HBK20 ‘Hubble bubble’ model on the same scales. This will 
help to assess whether the Hubble and bulk flow tensions faced by 
the ACDM model might have a common Milgromian solution. More 
generally, our results will help to clarify whether the velocity field 
in a local supervoid solution to the Hubble tension might also match 
the observed bulk flow on scales of up to 250 h~! Mpc. 

The structure of this paper is as follows: In Section 2, we describe 
how we obtain the predicted bulk flows in the same manner as the 
reported observations. We then present our results in Section 3 and 
conclude in Section 4. 


2 METHODS 


The starting point for our analysis is the peculiar velocity field of 
the semi-analytic void model of HBK20, which is shown in their 
figure 8 for the Maxwell—Boltzmann void density profile (velocity 
fields of the relevant models are shown in our Appendix A). This 
is a combination of a spherically symmetric outflow from the void 
centre with a systemic motion of the whole void towards the left, 
which reduces the spherical symmetry to axisymmetry. We refer to 
the peculiar velocity in the CMB frame as the total velocity Vot 
(analogous to equation 72 of HBK20). This is well suited for a 
comparison with observations because in studies of the bulk flow 
on large scales, heliocentric redshifts are typically corrected for the 
precisely known velocity of the Sun in the frame of the CMB.” 


2.1 Possible locations for the Local Group 


The position of the LG in this framework can be deduced by finding 
points which move with a velocity of viot = vig = 627 4 22 kms“! 
(Kogut et al. 1993), where we use the notation that v = |v| for any 
vector v. Fig. 8 of HBK20 shows the locus of such points with a solid 
black curve. Those authors suggested that the LG is towards the right 
hand side of this curve, which is located further away from the void 
centre. However, their study lacked any reliable way to precisely 
determine where the LG ought to be. 

We simplify our analysis by only considering two possible LG 
locations, which we call the inner and outer vantage point (relative 
to the void centre). These bracket the range of possible LG distances 
from the void centre. Importantly, both vantage points place the LG 
on the symmetry axis of the problem, allowing us to work with an ax- 
isymmetric code. We consider both vantage points separately for all 
three considered void (under)density profiles (Maxwell—Boltzmann, 
Gaussian, and Exponential). We also consider uncertainties in the 
model predictions for each of these six cases due to the uncertainty 
in upg, which can slightly vary the location of the LG by a few Mpc 
(see Appendix A). This is however less significant than observational 
uncertainties in the bulk flow. 

In Table 1, we show the possible locations of the LG for all three 
void density profiles. In all cases, the inner vantage point lies within 
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Table 1. The distance of the LG from the void centre (in Mpc) in each of the 
considered void density profiles for the two possible vantage points consistent 
with the CMB-frame velocity of the LG (see the text). The impact of its 22 km 
s7! uncertainty (Kogut et al. 1993) is considered separately in all cases. We 
also consider setting vug = 840 km s7! (Section 3.1), in which case the inner 
vantage point’s distance from the void centre is 115.99 (158.39) Mpc for the 


Gaussian (Exponential) void profile. 


Density profile 

LG location vig (km s7!) MB Gauss Exp 

Inner 627 + 22 136.38 142.59 192.19 
vantage 627 138.08 145.79 196.39 
point 627 — 22 139.78 149.29 200.69 
Outer 627 + 22 265.37 505.95 848.56 
vantage 627 262.27 494.76 816.16 
point 627 — 22 259.17 483.86 786.46 


200 Mpc of the void centre. The outer vantage point is generally 
much more distant. Only the Maxwell—Boltzmann profile gives an 
outer vantage point reasonably close to the void centre, at a distance 
of dup, = 262.27 Mpc. 


2.2 Observed peculiar velocity 


We describe the problem using Cartesian coordinates centred on 
the observer. Only two axes are needed because the problem is 
axisymmetric. The x-axis corresponds to the symmetry axis of 
the problem, while the y-axis lies in the orthogonal direction. The 
distance of any point from the observer is thus r = ,/x? + y?. 
Observations of distant galaxies can only tell us their peculiar 
velocity along the line of sight. Applying this consideration, the 
observable component of the peculiar velocity at point i is 


XUx + yvy 
Vobs,i = (==>) , 2) 


r 


where viot = (vx, vy) is the velocity in the CMB frame. This is found 
by adding the velocity relative to the void centre with the systemic 
velocity of the void as a whole, which is towards —x. The systemic 
void velocity is given in tables 4 and Cl of HBK20 and is not 
adjustable in our analysis. 


2.3 Bulk flow velocity 


While the line-of-sight peculiar velocity of point iis typically thought 
of as a scalar quantity, it will help to think of it as a vector Vops j 
pointing from the observer towards point i. The bulk flow velocity in 
the CMB frame is the weighted average of these line-of-sight peculiar 
velocity vectors within a spherical volume of radius py, centred on 
the observer. This definition matches equation 26 of Watkins et al. 
(2023). 

Due to the simulated velocity field being axisymmetric with 
respect to the x-axis and the observer also lying on this axis, the bulk 
flow must be along it. We therefore consider only the x-component 
of Vops,i; Which is given by 


Vobs,x,i = Vobs,i (=) . (3) 
PLE 


Calculating the bulk flow velocity is then just a matter of taking a 
suitably weighted average of this quantity, which we discuss below. 

Each point has a weight w; « V;/ r°, where V; is the volume covered 
by cell i. The factor of 1/7? is necessary because we are dealing with 
the radial velocity components instead of the 3D velocities. Under the 
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assumption that the velocity field is curl-free, this will then match the 
average of the 3D velocities, which is how the bulk flow is defined 
(see section 3 of Peery, Watkins & Feldman 2018). To minimize 
numerical issues at very low radii caused by the 1/7? factor, we 
soften it by adding to r in quadrature half the radial width of each 
cell used in our calculation. Due to the high resolution used, this has 
a negligible impact on our analysis. The bulk flow is then 


X; Wi Vobs,x,i 

Bw 
where the sums are taken over all points within 7pux of the observer. 
The volume weighting accounts for the inhomogeneous distribution 
of points, a consequence of the discretization scheme detailed in 
section 3.3.6 of HBK20 and reused here for simplicity. This scheme 
reduces the computational cost by working with a polar coordinate 
system centred on the void centre. Because we need to work with 
spherical regions centred on the LG rather than the void centre, we 
substantially increase the resolution compared to that used in HBK20. 
We also compute the velocity field out to a larger distance because 
we need predictions out to 250 h~! Mpc from the outer vantage point. 

Since the simulated velocity field is smooth, setting Fouk —> 0 
implies that all considered points have the same Vio. This velocity 
must be parallel to the symmetry axis as the velocity component 
perpendicular to it must become negligible at small distances from 
the axis. Equations (2) and (3) can then be combined and simplified 
to give 


(4) 


Vbulk (< Foulk) = 


xN? 2 
Vobsxi =V (7), = vecos a, (5) 
i 


where v = |Vial is the speed of any particle in the small region under 
consideration, œ; is the angle at the observer between the direction 
towards point 7 and the positive x-axis, and the power of 2 accounts 
for the projection of Vio onto the line of sight and the projection of 
the line of sight onto the direction of v,.,. Relating this to equation 
(4) shows that the observed bulk flow in this small region is only 
v/3 because the angle-averaged value of cos*a is 1/3 when there are 
three spatial directions. Note also that in a small sphere centred on 
the observer, we must have that Viot —> VLG aS Fouk —> O because 
we expect that the peculiar velocity of the LG arises mostly from 
structure on quite large scales. In the model, vyg is towards —x for 
the inner vantage point and towards +x for the outer vantage point. 
This is because the systemic velocity of the void dominates near its 
centre, while outflow from the void becomes more important further 
out. 


3 RESULTS AND DISCUSSION 


The bulk flow curves for all considered void profiles and LG locations 
(inner/outer) are shown in Fig. | for the nominal vyg. Results for the 
Maxwell—Boltzmann, Gaussian, and Exponential void profile are 
shown in green, red, and blue, respectively. Results for the inner 
(outer) vantage point are shown using a solid (dotted) line. Two 
additional curves are shown for the inner vantage point in which via 
is varied from its nominal value by its observational uncertainty: the 
dashed—dotted (dashed) curves assume a lower (higher) upg. The 
resulting uncertainty on the predicted bulk flow is much smaller than 
in the observed bulk flow amplitude (solid black points with error 
bars; Watkins et al. 2023). All the simulated bulk flow curves start at 
vic/3 for reasons discussed below equation (5). 

The bulk flow represents an average of the motions within a region. 
If the velocity field were linear in the position, then the bulk flow 
curves would remain flat at v_g/3. Non-linearity in the velocity field 
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causes the bulk flow curve to deviate from this value. To better 
understand the simulated bulk flow curves, it will help to consider 
the velocity field (Appendix A). 

We first consider the case where the LG is at the inner vantage 
point. The bulk flow curve for the Maxwell—Boltzmann profile is 
fairly flat because vio is nearly linear: if we consider its value along 
the symmetry axis for simplicity, it drops from +1500 km s! at the 
void centre to zero at ~200 Mpc (as indicated by the black star) and 
then rises back up again at a similar rate, albeit with a velocity in the 
opposite direction (Fig. Al). This causes the average motion within 
a sphere to stay roughly flat as we increase the radius of the sphere, 
leading to a fairly flat bulk flow curve. Since the observed bulk flows 
are much faster than vy g/3 and increase with scale, it is not possible 
for the Maxwell—Boltzmann profile to explain the observations (solid 
green curve in Fig. 1). 

The situation is rather different for the Gaussian and Exponential 
profiles, which yield bulk flow curves that rise rather quickly. If we 
again consider Vo along the symmetry axis, it rapidly drops from 
~2000 km s7! at the void centre to 627 km s7! at the location of the 
LG, but v; then only very gradually decreases to zero (see Figs A2 
and A3, respectively). This causes the average motion within a sphere 
to rise quickly with the radius of the sphere, leading to a rising bulk 
flow curve. We expect the bulk flow to remain parallel to vig at all 
radii, which is approximately the case (Kogut et al. 1993; Watkins 
et al. 2023). 

We find that both the Gaussian and Exponential density profiles re- 
sult in bulk flow curves with a similar trend to the observed bulk flow 
curve: there is an initial steep increase which then almost stagnates 
at large radii. However, both models undershoot the observations by 
50-100 km s~!, with the largest deviation at ry œ~ 150h7! Mpc. 
Part of this can be attributed to relativistic corrections, which imply 
that the actual bulk flows are about 10 percent smaller than reported 
(Heinesen 2023). However, these corrections scale with z and are in 
any case too small to fully explain the discrepancy. 


3.1 Our vantage point 


The bulk flow around the outer vantage point completely differs from 
the inner one, as shown by the continuously declining thin dotted 
curves in Fig. 1 (the same colour is used for each density profile). 
If the LG were located at the outer vantage point, then regardless of 
the void density profile, the bulk flow would decrease with rpux 
and reverse direction at buk © 30047! Mpc (420h7! Mpc) for 
the Maxwell—Boltzmann (Gaussian) density profile. The bulk flow 
curves for the outer vantage point are in catastrophic disagreement 
with the observed fact that Vpux >> ViG/3 for Fouk œ% 200 A7! Mpc. 
The outer vantage points are at least 2x further away from the void 
centre than their respective inner counterparts (4x for the Gaussian 
and Exponential profiles). Thus, our results seem to imply that we 
need to be within ~200 Mpc of the void centre (see Table 1), in 
agreement with the conclusion of HBK20. Being fairly close to the 
centre of a void would also limit the extent to which the Universe 
appears anisotropic. 

The impact of the 22 km s~! uncertainty in vrc is illustrated in Fig. 
1, where the dashed (dashed—dotted) lines show the impact of raising 
(reducing) upg by this amount from its nominal value of 627 km svt. 
The uncertainty in vig has little impact on the observed bulk flow 
for the inner observers because the modified velocity of the LG only 
slightly shifts its position (see Table 1 and Appendix A). The resulting 
changes to Vpux are <10 km E Increasing vrc by its uncertainty 
causes the bulk flow curves for the Gaussian and Exponential profiles 
to slightly shift towards the observational data, while decreasing 
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Figure 1. Bulk flows in spheres of different radii around an observer with a suitable CMB-frame velocity for the Maxwell—Boltzmann, Gaussian, and Exponential 
density profiles (green, red, and blue curves, respectively). The solid black points with uncertainties show the observed bulk flows (bottom right panel of figure 7 
in Watkins et al. 2023). The two innermost data points use the same method and survey data but were not published in that study, so these were provided in a 
private communication by its lead author. At roux = 0, the simulated curves begin at v_G/3 because the observed bulk flows use only the line-of-sight velocities, 
which are treated as vectors that are then combined (see below equation 5). The solid curves represent the bulk flows as seen by the inner observer if vg = 
627 km s~!. The dashed curves (typically above each solid curve) are the bulk flows for vig = 627 + 22kms~!. The dashed—dotted curves (typically below 
each solid curve) are the bulk flows for vua = 627 — 22 km s™!. These additional curves help to show the impact of the uncertainty in vig, which slightly shifts 
our vantage point in the context of each model and can thus affect the predicted bulk flow. The grey dashed (dotted) curve shows the bulk flow curve for the 
Gaussian (Exponential) density profile if vug = 840 km s~!, which shifts the inner vantage point by 29.8 (38.0) Mpc towards the void centre. The coloured thin 
dotted curves represent the outer observers for all density profiles. These all assume the observed vig and use the same colour as the corresponding results for 


the inner observer. 


vic has the opposite effect. However, these differences are barely 
perceptible. 

The bulk flow curves from the Gaussian and Exponential density 
profiles underestimate the observed bulk flow curve. Note that since 
the bulk flow measurements are already in the CMB frame, the 
assumed vyg enters our analysis only by setting our location within 
the void. The closer the vantage point is to the centre of the void, the 
higher the bulk flow curve becomes overall. Considering the inner 
vantage point with any density profile, the uncertainty in vyg of 
22 km s~! corresponds to a shift to our position by a mere 3—4 Mpc. 
As discussed earlier, this is barely perceptible. 

If we allow a larger positional shift of the LG, we could in 
principle match the observed bulk flow curve fairly well. Such an 
increase to upg is plausible because ideally our location within 
the void would be set using an estimate of Vioca, the bulk flow 
on a scale large enough to average over motions at the scale of 
individual galaxy groups and clusters but small compared to the 
KBC void. We assumed vjocaj = Vig as an approximation due to the 
LG frame being the largest reference frame to which a transformation 
can easily be performed. However, the true value of vica will 
generally differ from upg. The investigated model by design ignores 
the possibility of small-scale structures as it only considers the 


MNRAS 527, 4388—4396 (2024) 


much larger structure that is the KBC void. In the real Universe, 
additional structures like gravitational attractors exist in the form 
of galaxy clusters, for instance the Virgo and Fornax clusters at 
distances of 17 and 20 Mpc, respectively (Mei et al. 2007; Blakeslee 
et al. 2009). Small-scale underdensities can also play a role. In the 
local vicinity of such small-scale structures, the velocity field would 
become distorted, thereby shifting the appropriate vantage point. 
The larger scale velocity field of the void would remain unaffected 
by this. 

Interestingly, the bulk flow on a scale of 50h~! Mpc is observed 
to be 297 + 14 km s™! (private communication from R. Watkins; 
same method and data as in Watkins et al. 2023). This suggests an 
actual bulk motion on this scale of 891 + 42 km s~! based on our 
argument in equation (5). Allowing for the fact that the predicted 
bulk flow curve rises by just under 15 kms~! between ru = 0 and 
Foulk = 50 A7! Mpc for the most realistic void models from HBK20, 
the above result suggests that an appropriate choice of Vioca is slightly 
below 891 km s~!. This motivates us to consider the case vig = 
840 km s~! and update our vantage point accordingly. The results are 
shown using the dashed (dotted) grey line in Fig. 1 for the Gaussian 
(Exponential) profile. In both cases, the model agrees quite well 
with the observed bulk flow curve, perhaps indicating that the local 
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velocity due to large-scale structure is more accurately reflected by 
the bulk flow on a 50/7! Mpc scale than by vic. Regardless of the 
considered void models, the mismatch between the observed vig 
and 3vpyx suggests that the LG is moving +200 km s~! slower in 
the CMB frame than is typical for the matter within 507! Mpc. 
This modest difference could plausibly arise from nearby structures 
beyond the LG. 


3.2 Predictions at large radii 


For all considered density profiles, the HBK20 model predicts a 
rising bulk flow amplitude over the range rpux = 250 — 500 A7! Mpc 
if we are located at the inner vantage point, as strongly suggested 
by our results. At 50047! Mpc, the Gaussian and Exponential 
profiles reach bulk flow amplitudes of approximately 475 and 520 km 
s—!, respectively. However, these predictions should be treated with 
caution because it is likely that our simplified model for the KBC void 
breaks down at such large distances as there must be other structures 
in the Universe. Even within the model we are using, it is necessary to 
assume more distant structures that impose an external gravitational 
field on the void as a whole (HBK20). At some point, it is no longer 
possible to approximate the structures responsible for this external 
field as being much more distant than the region under consideration. 
These issues should be more thoroughly addressed using a Gpc-scale 
cosmological simulation of the vHDM model and analysis of KBC 
void analogues therein (Russell et al., in preparation). It will be 
interesting to see if the bulk flow typically reverses on larger scales 
than those considered by Watkins et al. (2023), as would be needed to 
explain the dipole anisotropy in the galaxy cluster distance-redshift 
relation (Migkas et al. 2021). 

A local void solution to the Hubble tension necessarily implies 
that the Planck cosmology works at suitably high redshift. While 
our focus here is on the predictions of the HBK20 model at the 
present epoch, those authors used the lightcone analysis described 
in their section 3.3.3 to show that the void-induced enhancement to 
the redshift decays rather slowly and is not negligible 8 Gyr ago 
(see their figure 16). This corresponds to z © 1, which can seem 
far beyond the reach of a local void. But if the void is fairly large 
and given also the rather slower decline of MOND gravity with 
distance compared to Newtonian gravity, it is quite possible that the 
impact of the void on cosmological observables is not restricted to 
z < 0.1, as was assumed in figure 1 of Kenworthy, Scolnic & Riess 
(2019). It is thus not surprising that those authors were unable to 
find any evidence of a local void (other problems with a local void 
scenario were addressed in detail in section 5.3 of HBK20). It is 
interesting in this regard that the data of Kenworthy, Scolnic & Riess 
(2019) only goes up to z = 0.5 (see their figure 5). More recent 
studies do find a change to the inferred Ho at higher redshift if the 
data are binned in z and care is taken to infer the cosmological 
parameters from only the data within each redshift bin, removing 
residual correlations with the data in other bins (Jia, Hu & Wang 
2023). Their analysis clearly shows a return to a Planck cosmology at 
higher redshift, regardless of whether we consider supernovae alone 
or combine it with other data sets to improve the accuracy. A similar 
result is apparent in the preliminary analysis of Gémez-Valent et al. 
(2023), whose reconstructed expansion rate history closely follows a 
Planck cosmology at high redshift and deviates only when z < 0.6. 
We are currently exploring the consistency of the HBK20 model 
with cosmological data sets at high redshift, building on its known 
success fitting supernova data out to z = 0.15 (Mazurenko et al., in 
preparation). 
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4 CONCLUSIONS 


It has been suggested that the anomalously high local Ho could 
be due to outflows from a local void (HBK20). We revisit their 
best-fitting semi-analytic void model for the three initial (z = 9) 
density profiles they considered (Maxwell—Boltzmann, Gaussian, 
and Exponential). In each case, we consider two possible vantage 
points where the simulated peculiar velocity in the CMB frame 
matches that of the LG. For these six combinations of void density 
profile and LG location, we determine the bulk flow vpux in spheres 
of different radii reux < 250 A7! Mpc. This predicted bulk flow curve 
is compared with recent observations (Watkins et al. 2023) based on 
the CosmicFlows-4 galaxy catalogue (Tully et al. 2023), which is 
designed to obtain reliable distances and thus line-of-sight peculiar 
velocities. We also consider the small uncertainty in the predicted 
bulk flow due to the uncertainty in vig, which slightly shifts the 
location of the LG. The impact on the predicted vp, is much smaller 
than its observational uncertainty (Fig. 1). 

With all three density profiles, the predicted bulk flow curve 
from the outer vantage point (relative to the void centre) is strongly 
incompatible with observations, both qualitatively and quantitatively. 
Assuming instead that we are located at the inner vantage point, the 
Maxwell—Boltzmann density profile also gives a bulk flow curve that 
is observationally excluded (albeit in somewhat better agreement). 
However, the Gaussian and Exponential profiles give bulk flow 
curves that qualitatively agree quite well with the observational 
data, nicely capturing the steep increase followed by a flattening 
out. Within the context of the HBK20 model, this implies a void 
which is deepest in the centre and that we are located within 200 
Mpc of the void centre while remaining consistent with the observed 
LG peculiar velocity. Taking into account its 22 km s~! uncertainty 
barely impacts the bulk flow curve. However, if one allows for small- 
scale flows within the void which are not captured by the model 
(as discussed in Section 3.1), then a good quantitative agreement 
is achievable (grey curves in Fig. 1). These curves nicely match 
the bulk flow measurement of Watkins et al. (2023) on a scale 
of 50h47! Mpc, suggesting that the velocity of the LG has been 
reduced slightly by local structures on a scale not considered in the 
model. 

The agreement outlined above is reasonable considering the 
simplicity of the semi-analytic model used, which assumes spherical 
symmetry relative to the void centre plus a constant extra velocity 
everywhere within the void due to the impact of even more distant 
structures (HBK20). Obviously the local Universe does not have 
such a simplified velocity field, which must instead be tied to its 
complex history of structure formation over a Hubble time. This 
should be explored using a Gpc-scale cosmological simulation of 
the vHDM paradigm, building on previously published smaller scale 
simulations (e.g. Katz et al. 2013; Wittenburg et al. 2023) — but 
ideally reaching scales similar to the MXXL simulation in the ACDM 
framework (Angulo et al. 2012). This will address the plausibility 
of the observed bulk flows in a self-consistent vHDM cosmological 
simulation (Russell et al., in preparation). In addition, detailed future 
observations of the bulk flow beyond 250 h7! Mpc can be compared 
to the bulk flows predicted by the here utilized models and those 
found in Gpc-scale vHDM simulations. Once such measurements 
become available, the tension with the ACDM model should also be 
quantified. It is quite possible that future data will worsen the 4.80 
tension reported by Watkins et al. (2023) on a scale of 200h—! Mpc, 
especially given that the bulk flow seems to be increasing with scale 
while the typically expected value in ACDM decreases with scale 
(see their figure 7). 
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It is important to stress that the recent bulk flow measurements 
of Watkins et al. (2023) were not available to HBK20, who did not 
in any way constrain the void parameters using the observed bulk 
flow. Despite this, their model naturally accounts for the observed 
rising behaviour of the bulk flow on scales of 100 — 250 A7! Mpc 
without any parameter adjustments for two of the six considered 
combinations of void profile and observer location, which was chosen 
to match the peculiar velocity of the LG. This agreement suggests 
that a Milgromian cosmological framework can consistently explain 
the Hubble and bulk flow tensions by assuming that we are located 
in a large supervoid similar to the observed KBC void. Our analysis 
also imposes constraints on the position of the LG relative to the 
KBC void centre, and to a lesser extent on the void density profile. 
Future work will test this model through better measurements of 
the KBC void density profile and the bulk flows on larger scales 
than those probed so far. In addition, it will be important to explore 
self-consistent cosmological simulations rather than the simplified 
semi-analytic model considered here. 

Our results suggest that many of the severe cosmological tensions 
faced by ACDM could be alleviated if structure grows more rapidly 
than it predicts. This could also explain the unexpectedly large 
integrated Sachs—Wolfe effect associated with supervoids (DES 
Collaboration 2019; Kovacs et al. 2022) and the 2 percent dipole 
in the distribution of quasars (Secrest et al. 2021, 2022; Dam, 
Lewis & Brewer 2023). In this scenario, overdensities should also 
be more pronounced than expected in standard cosmology. There is 
strong evidence for this from massive high-redshift galaxy cluster 
collisions like El Gordo, which has been shown to falsify ACDM at 
>5o confidence for any plausible infall velocity (Asencio, Banik & 
Kroupa 2021, 2023). Even the formation of individual galaxies could 
be faster than expected in ACDM, as suggested by several recent dis- 
coveries with the JWST of massive galaxies in the very early Universe 
(z 2 10; Haslbauer et al. 2022, and references therein). Such rapid 
galaxy formation can be naturally understood in MOND (Eappen 
et al. 2022). While it is too early to say what might be responsible for 
these tensions with ACDM, the wide range of different observational 
techniques used to reveal these inconsistencies makes it likely that 
many if not all of them are genuine problems for the paradigm that 
go beyond the background expansion history. 
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APPENDIX A: VOID VELOCITY FIELDS 


While our focus has been on the bulk flow, it can be helpful to consider 
the velocity field in the void models we explore. We therefore 
show a 2D map of Vio for the Maxwell—Boltzmann, Gaussian, and 
Exponential density profiles in Figs A1,A2, and A3, respectively. 
The velocity fields are axisymmetric around the x-axis, with the 
void as a whole moving towards the left. The velocity fields are 
shown in polar coordinates centred on the void centre, which is 
situated at ryoiq = 0. While the direction of viot is not shown for 
simplicity, this is directed approximately radially away from the 
point where vio. = 0, which we indicate with a black star. In all these 
figures, the solid black line shows the locus of points where vi. = 
vto = 627 km s~!, while the dashed and dashed—dotted black lines 
represent the curves along which viot = 627 + 22 km s7! and Vot = 
627 — 22 km s7!, respectively, thus bracketing the observational 
uncertainty in upg. The grey line represents vyg = 840 km s“!, 
which is motivated by the bulk flow measurement of Watkins et al. 
(2023) on a scale of 50 h—! Mpc (see Section 3.1). The lower panels 
show a close-up of the region where viot = 627 + 22 km s~!. The 
vantage points considered for the bulk flow calculations are at the 
two intersections of each isovelocity contour with the x-axis. Our 
results strongly suggest that the LG is within 200 Mpc of the void 
centre. 
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Figure Al. Map of the total velocity viot in the CMB frame resulting from 
the Maxwell—Boltzmann void density profile. Results are axisymmetric about 
the x-axis because the model applies a Galilean transformation to a spherically 
symmetric outflow (HBK20). Distances and angles shown here are relative to 
the void centre. The solid black curve shows the locus of points where vtot = 
utc. To highlight its uncertainty, the dashed—dotted (dashed) line shows the 
locus of points where vig is assumed to be smaller (larger) than its observed 
value by its observational uncertainty. The grey line corresponds to vig = 
840 km s7! (Section 3.1). The bottom panel shows a zoomed in view. Vtot is 
directed roughly away from the black star in both panels marking the point 
where vtot = 0. 
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Figure A2. Similar to Fig. Al, but for the Gaussian density profile. Notice Figure A3. Similar to Fig. Al, but for the Exponential density profile. 
that the isovelocity contours are elongated towards the right, indicating that 

the velocity field changes more slowly at larger distances from the void centre. This paper has been typeset from a TEX/IATEX file prepared by the author. 


This causes the rising bulk flow curve evident in Fig. 1 (see Section 3). 
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